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We rewrite the complex Klein-Gordon (KG) equation with a mexican-hat scalar field potential 
in a thermal bath with one loop contribution as a new Gross-Pitaevskii-like equation. We interpret 
it as a charged and finite temperature generalization of the Gross-Pitaevskii equation. We find its 
hydrodynamic version as well and using it, we derive the corresponding thermodynamics. We obtain 
a generalized first law for a charged Bose-Einstein Condensate (BEC). We translate the breaking of 
the U(l) local symmetry of the KG field into the new version of the Gross-Pitaevskii equation and 
demonstrate that this symmetry breaking corresponds to a phase transition of the gas into a BEC 
and show the conditions for condensation and/or phase transition for which this system naturally 
becomes superfluid and/or superconductor. 
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INTRODUCTION 

Since its observation with the help of magnetic traps 
[H, the phenomenon of Bose-Einstein condensation has 
spurred an enormous amount of works on the theoretical 
and experimental realms associated to this topic. The 
principal interest in the study on Bose-Einstein conden- 
sation is its interdisciplinary nature. From the thermo- 
dynamic point of view, this phenomenon can be inter- 
preted as a phase transition, and from the quantum me- 
chanical point of view as a matter wave coherence aris- 
ing from overlapping de Broglie waves of the atoms, in 
which many of them condense to the grown state of the 
system. In quantum field theory, this phenomenon is re- 
lated to the spontaneous breaking of a gauge symmetry 

. Symmetry breaking is one of the most essential con- 
cepts in particle theory and has been extensively used 
in the study of the behavior of particle interactions in 
many theories Q . Phase transitions are changes of state, 
related with changes of symmetries of the system. The 
analysis of Symmetry breaking mechanisms have turn out 
to be very helpful in the study of phenomena associated 
to phase transitions in almost all areas of physics. Bose- 
Einstein Condensation is one topic of interest that uses 
in an extensive way the Symmetry breaking mechanisms 

01 , its phase transition associated with the condensation 
of atoms in the state of lowest energy and is the con- 
sequence of quantum, statistical and thermodynamical 
effects. 

On the other hand, the results from finite tempera- 
ture quantum field theory 0, [B[ raise important chal- 
lenges about their possible manifestation in condensed 
matter systems. By investigating the massive Klein- 
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Gordon equation (KG), in [6fl we were able to show, that 
the KG equation can simulate a condensed matter sys- 
tem. In 6] it was shown that the KG equation with a self 
interacting scalar field (SF) in a thermal bath reduces to 
the Gross-Pitaevskii equation in the no-relativistic limit, 
provided that the temperature of the thermal bath is 
zero. Thus, the KG equation reduces to a generalized 
relativistic, Gross-Pitaevskii (GP) equation at finite tem- 
perature. But a question remains open. The KG equa- 
tion with a self interacting SF potential defines a sym- 
metry breaking temperature, at which the system exper- 
iments a phase transition. However, this phase transition 
does not necessarily means a condensation of the parti- 
cles of the system. 

In the present work we study the complex KG equa- 
tion with a Mexican hat SF potential in a thermal bath. 
The idea is very simple, the KG equation up to one loop 
in perturbations is able to explain the phase transition 
of a SF, like the Higgs field, close to the moment of the 
phase transition, when the KG equation breaks the U(l) 
symmetry of the corresponding Lagrangian. On the other 
hand, the Gross-Pitaevskii equation is able to explain the 
behavior of a BEC at zero temperature. In this work we 
wonder if the KG equation in a thermal bath is able to 
generalize the Gross-Pitaevskii equation in a region close 
to the phase transition of the Bose gas into a BEC, be- 
cause the KG equation contains the information of the 
temperature of the bath in the SF potential. In base of 
this idea, we rewrite the complex KG equation into a 
Gross-Pitaevskii-like equation an interpret it as a Gross- 
Pitaevskii one at finite temperature. In order to see if 
this equation can explain the phase transition of a Bose 
gas into a BEC we write its corresponding thermodynam- 
ics, deriving a corresponding first law for the Bose gases. 
The only difference we find here with respect to the tradi- 
tional first law of the thermodynamics is a therm where 
the quantum mechanical character of the KG equation 
is present. On the other hand, a phase transition does 
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no necessarily means gas condensation. Therefore we ob- 
tain the temperature and the conditions of condensation 
of this Bose gas. At the end of the work we qualita- 
tively see under which conditions the Bose gas becomes 
supcrfluid or/and superconductor. 



GAUGE SYMMETRY BREAKING 

We start with a model having a local U(l) symmetry 
given by the lagrangian 

C = - (V„$* + ieA^*) (V"$ - ieA»$) 



(i) 



where F M „ = V^-A^ — V ' is the Maxwell tensor and 
the SF potential V is the easiest case of a double- well in- 
teracting Mexican-hat potential for a complex SF $(a;, t), 
interacting inside a thermal bath in a reservoir that can 
have interaction with its surroundings up to one loop of 
correction, that goes as 
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where m 2 — m 2 c 2 /K 2 is the scalar particle mass, A = 
X/(h 2 c 2 ) is the parameter describing the interaction, k B 
is the Boltzmann's constant, h is the Planck's constant, c 
is the speed of light and T is the temperature of the ther- 
mal bath, this result includes both quantum and thermal 
contributions. 

From here the interpretation of the different quanti- 
ties of the system is clear, $ is the KG version of the 
Gross-Pitaevskii function ^ (see bellow), the Maxwell 
field is an electromagnetic field in the system, A is the self- 
interaction parameter related with the scattering length 
of the system and T is the temperature of the thermal 
bath where the system lie. 

The dynamics of a SF is governed by the KG equation, 
it is the equation of motion of a field composed of spinless 
particles. In order to confine the BEC experimentally 
one needs to add an external field controlled by hand in 
order to cause the condensation, like a laser for example. 
In order to do this, we will add an external field <fr that 
will interact with the SF to first order, such that the KG 
equation will be given by 



9 dV , 
□ p$ 2to 2 <M> = 0, 



(3) 



where A^ is the electromagnetic four potential. Observe 
that we can rewrite the D'Alambertian as 

□| = (V-2kA)-V-(~-2iep)~-ieV M A*-e 2 V4'« 
c at cot 

(5) 

where we have written the electromagnetic four potential 
as A M = (A, ip). In what follows we will use the Lorentz 
gauge V M A M = 0. It is convenient to consider the total 
(effective) potential V adding the external one and the 
term e 2 A^A^ = e 2 A 2 to the potential, such that 
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where we see that the parameter e is the coupling con- 
stant between the electromagnetic and the SFs. In terms 
of the effective potential Vr the KG equation now can be 
written as 



dVr 



0. 



(7) 



where now O 2 = V M V M . The Maxwell equations also 
read 

V M F — — i 

= ie(<rV I ,$-$V I ,$*) + 2e 2 <M>*v4 i ,. (8) 

We can define an effective mass by 

m e f f c 2 = y/m 2 c 4 + e 2 A 2 . (9) 

where e = ehc. For the Vr potential ©, the critical 
temperature T c s B where the minimum of the potential 
$ = becomes a maximum and at which the symmetry 
is broken is, 

2c 2 / 
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Potential (J6j> has a minimum in $ = when the temper- 
ature T > Tf B . If T < T^ B , the point $ = becomes a 
maximum and potential (JBJ) has two minima in 
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being $ = Re t& . In the maximum $ = the second 
derivative of the potential Vr with respect to the SF reads 
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where for a charged field the D'Alambertian operator is 
given by, 



□ | = (V M - ieA M ) (V - ieA") 



(4) 
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In what follows we will re-write the KG equations in 
order to interpret the KG equation as a Gross-Pitaevskii 
one. 



THE GENERALIZED GROSS-PITAEVSKII 
EQUATION 

Now for the SF we perform the transformation 

In terms of the complex function , the KG equation 
d2]) now reads, 
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where we have kept just the equation for the \t part, 
the complex conjugate can be described in the same way. 
The notation used is ' = 1/cd/dt (fT3")l is the KG equation 
(|3]) or ([7]) rewritten in terms of the function ^ and tem- 
perature T. This equation is an exact equation defining 
the field ^(x, t), where <\> defines the external potential 
acting on the system and the term with A represents the 
self-interaction of the SF within the system. We will 
consider equation (fT3|) as a generalization of the Gross- 
Pitaevskii equation for finite temperatures and relativis- 
tic particles. This is because when T = and in the 
non-relativistic limit, D 2 — ► V 2 , eq. (IT31 becomes the 
Gross-Pitaevskii equation for Bose- Einstein Condensates 
(BEC), provided that A = 8nh 2 c 2 a k 2 , being a the s-wave 
scattering length @ . The static limit of equation (IT3l) is 
known as the Ginzburg-Landau equation. Observe that 
the temperature T and the external field <fr must be ma- 
nipulated from outside. 

At this point it is noteworthy to mention the order 
of magnitude of the quantity pop, the associated crit- 
ical temperature at which the symmetry of the sys- 
tem is broken. Assuming that the external potential 
<j> and the additional field A 2 are time independent, 
i.e, that the system is in static thermal equilibrium, 
we can write ksT^ B ~ in the center of the sys- 

tem. If we set A = 8nh 2 c 2 an 2 into (10]), thus T c ss ~ 

Joules- 1 ™-*/ 2 ("A /cm x 10 62 ^ p or mstance m 
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the case of 87 Rb, with a mass m 
and a scattering length a ~ 5.2 x 10~ 7 m, the crit 
ical temperature of symmetry breakdown of this sys- 
tem is T c ss ~ ^ x 10 44 K, which depending on the 
value of k, it could be very big. However, for ex- 
ample, if we set n ~ 5 x 10 50 Joules -1 m~ 3 / 2 (see 
also Q), we obtain that T$ B - 3 x K. Nev- 
ertheless, the density n does not depend on the value 
of k, because n = k 2 k B (T c ss ) 2 /4(1 - (T/T C SB ) 2 ) ~ 



K 2 k B (T c ss ) 2 /4 = m 2 c 2 /(8Trh 2 a), and this quantity for 
the 87 Rb is ~ 10 36 /cm 3 . Multiplying this density times 
the mass, we obtain the corresponding mass density, 
which is ~ 1.9 x 10 14 gr/cm 3 . This mass density cor- 
responds to the one of a nuclear atom or a neutron 
star. Furthermore, for this values the self-interaction 
parameter A is very big A ~ 3.2 x 10 43 . Let us sup- 
pose for a moment that neutrons can build cooper pairs 
in a neutron star. In that case, the mass of two neu- 
trons is m = 2 x 1.6 x 10 _24 gr, the scattering length 
of the neutrons is a ~ 10 cm. Thus, the critical 
temperature of symmetry breakdown of this system is 
t sb ^ MLiq^x. If we set the typical temperature of 
collapse for a neutron star T ~ 10 12 K as the temperature 
of the symmetry breakdown, we obtain that k ~ 10 32 
Joules -1 m~ 3 ' 2 . In this case A ~ 25, the particle density 
n ~ 3.2 x 10 37 /cm 3 and the mass density ~ 10 14 gr/cm 3 , 
which is again of the order of magnitude of neutron stars 
(see [3 ]■ Other interesting example is BECs as Dark 
Matter llj] . In this case the mass of the SF could be of 
the order of the axion mass m ~ 0.1eV= 1.78 x 10 _34 gr. 
Here there are two cases, if the SF is an axion [ijj], the 
self-interaction parameter is A ~ 10 15 [13|, for this value 
of the mass and the self-interacting parameter, the criti- 
cal temperature is T C SB - 7 x 10~ 5 K= 6 x 10" 9 eV. The 
mass density is as expected p ~ 10 _45 gr/cm 3 . On the 
other hand, if the SF are Scalar Field Dark Matter par- 
ticles [l3|, the mass is again m ~ 0.1eV= 1.78 x 10~ 34 gr, 
but the self- interaction parameter is A ~ 10~ 6 . For that 
the critical temperature is T C 5B - 2.3 x 10 6 K= 200eV for 
a mass density of the order of the critical density of the 
universe p ~ 10 _29 gr/cm 3 . 



THE HYDRODYNAMICAL VERSION 

In what follows we transform the generalized Gross- 
Pitaevskii equation (|13[) into its analogous hydrodynam- 
ical version, (IR Il6l | , for this purpose the ensemble wave 
function ^ will be represented in terms of a modulus n 
and a phase S as, 



(14) 



where the phase S(x,t) is taken as a real function. As 
usual this phase will define the velocity. Here we will in- 
terpret n(x,t) as the number density of particles in the 
symmetry broken state, such that k 2 |^| 2 = k 2 ^^* = n, 
where k is the scale of the system, which is to be de- 
termined by an experiment, being both, S and n, func- 
tions of time and position. So, from this interpretation 
we have that when the KG equation oscillates around 
the $ = minimum, the number of particles in the 
ground state is zero, n = = p. Below the crit- 
ical temperature T c , close to the second minimum, 
R 2 mm = /c|((T c SB ) 2 - T 2 )/4, the density will oscillate 



4 



around n — n 2 k 2 B ((T^ B ) 2 — T 2 )/4 as can be seen by equa- 
tion (|14[) . In order to see this more clear, we perform the 
Madelung transformation (fl4l in the generalized Gross- 
Pitaevskii equation (|13l) . 

From (I13[) and (| 14[) . separating real and imaginary 
parts we obtain 
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Taking the gradient of (|15b[) and using the definitions of 
the fluxes 
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equations (11 5|) can be rewritten as, 
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(16b) 
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(17) 

j = 0, (18a) 
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where E = —1/cdA/dt + V • <p and £? = V x A re- 
spectively are the electric and the magnetic field vec- 



tors. Notice that in (|18bl) h enters on the right-hand side 
through the term containing the gradient of n. This term 
is usually called the 'quantum pressure' and is a direct 
consequence of the Heisenberg uncertainty principle, it 
reveals the importance of quantum effects in interacting 
gases. Multiplying by n, (|18b[) can be written as: 

nv + n(v ■ V)v = nF E + nF^ -Vp + nF Q + Vct,(19) 

where Fe = z\ (cE + v x B) is the electromagnetic 
force, F<f, — — V0 is the force associated to the exter- 
nal potential 4>, p can be seen as the pressure of the 



SF gas that satisfies the equation of state p = wn 2 , be- 
ing uj — A/(2to 2 c 2 k 2 ) an interaction parameter. Vp are 
forces due to the gradients of pressure, Fq = — VE/q is 
the quantum force associated to the quantum potential, 
10, 
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m z 



(20) 



and Vcr is defined as 
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where the coefficient C is given by 
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Observe that using (|18ap the term V(lnn) can be written 
as 



hn 2 ( 1 • 
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System (|18p is the hydrodynamical representation to 
equation (|13D and up to one constant, it is completely 
equivalent to (H"3l) . 



THE NEWTONIAN (NON-RELATIVISTIC) 
LIMIT 



Neglecting second order time derivatives and products 
of time derivatives we can simplify system (1181) . In this 
limit we arrive to the non-relativistic system of equations 



h + V • (nv) = 0, (22a) 
nv + n(v-V)v = nF E + nF ^ — Vp + nFQ + Vcr. 

(22b) 



Equation (122ap is the continuity equation, and (|22bp is 
the equation for the momentum. Observe that this last 
one contains forces due to the external potential, to the 
gradient of the pressure, viscous forces due to the inter- 
actions of the condensate and forces due to the quantum 
nature of the equations. Quantity V(lnn) plays a very 
important roll, in this limit it reads 



V(lnn) = -V(V • v) - V[V(lnrc) • v}. 



Thus 
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where now we have 



c = - 



2m 



r V • (nv), 



We interpret the function Vtr as the viscosity of the sys- 
tem, it contains terms which are gradients of the tem- 
perature and of the divergence of the velocity and den- 
sity (dissipative contributions) . The measurement of the 
temperature dependence in this thermodynamical quan- 
tity at the phase transitions might reveal important in- 
formation about the behavior of the gas due to particle 
interaction. 



THE THERMODYNAMICS 

In what follows we will derive the thermodynamical 
equations for the non-relativistic limit from the hydro- 
dynamical representation. We can derive a conservation 
equation for a function a, starting with the relationship 

(na)' = na + an (24) 

where a can take the values of 4>, Uq, etc., all of them 



fulfil equation ([24| . Using the continuity equation (|22a[) 
in (|2"4"|) we obtain, 

(na)' + V • (nva) — —nv ■ F a + na. 

Nevertheless, this procedure is not possible for a because 
in general we do not know it explicitly, only in some cases 
it might be possible to integrate it. 

Observe how the quantum potential Uq also fulfills the 
following relationship 



iUq + V ■ (nv p ) = 0, 



(25) 



which follows by direct calculation, and where we have 
defined the velocity density v p by 

(Vlnn)', 
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which can be interpreted as a velocity flux due to the po- 
tential Uq. Using the continuity equation (|22a[) . equation 
?5l) can be rewritten as 



(nU Q ) + V • (nU Q v + J p ) +nv-F Q = (26) 
where we have defined the quantum density flux 



J p = nv p . 



Equation (|26p is another expression for the continuity 
equation of the quantum potential Uq. 

As we know, for the non-relativistic limit the total en- 
ergy density of the system e is the sum of the kinetic, 



potential and internal energies 18J, in this case we have 
an extra term Uq due to the quantum potential, 



e = 2 nv2 + n § + nu 71 Uq + V'e (27) 
being u the inner energy of the system and 



1pE = —{<P -v ■ A) 



(28) 



the electromagnetic energy potential, defined in terms 
of the vector potential A and the electric potential (p. 
Observe that ipE fulfills the continuity equation 

(nip E ) + V • (nvipE +3b)= nips - nv ■ F E (29) 

being j B given by the continuity equation of the vector 
potential A 

— + (v-V)A = -(A ■ V)v + -j B , (30) 



Then from (j27|) we have that u will satisfy the equation 

(nu) + V • J u - V • J p + n4> = -pV • v, (31) 

being J u the energy current, given by a energy flux and 
a heat flux, J q , 

J u = nuv + J q + J b — pv, 

where V • J q — v ■ (Vcr), and V ■ J b — v ■ (nj B ), expres- 
sions that as we can see are related in a direct way to the 
velocity and gradients of temperature in the condensate, 
and are the ones that show in an explicit way the tem- 
perature dependence of the thermodynamical equations. 
With these definitions at hand we have, 

(nu)'+^ -(nvu+J q +J B—pv—J p )+n(j) = —pV-v. (32) 

In order to find the thermodynamical quantities of the 
system in equilibrium (taking p as constant on a given 
volume), we restrict the system to the regime where the 
auto-interacting potential is constant in time, with this 
conditions at hand for (|3"2"|) we have: 

(nu) ' + V • (nvu + J q + J b — pv — J p ) = —pV ■ v (33) 

From (|33[) we can have a straightforward interpretation 
of the terms involved in the phase transition. As always 
the first term will represent the change in the internal 
energy of the system, — pV • v is the work done by the 
pressure and V • v is related to the change in the volume, 
J q contains terms related to the heat generated by gradi- 
ents of the temperature VT and dissipative forces due to 
viscous forces ~ V(V-i>) and finally but most important 
we have an extra term, V ■ J p , due to gradients of the 
quantum potential (|20[) . 
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Integrating this resulting expression on a close region, 
we obtain 

J nudV + <j> \j q + Jb + pv) ■ ndS - <j> J p -ndS 



- P dt I ^ 



Equation (|33[) is the continuity equation for the internal 
energy of the system and as usual, from here we have 
an expression that would describe the thermodynamics 
of the system in an analogous way as does the first law 
of thermodynamics, in this case for the KG equation or 
a BEC. This reads 



dU = AQ + dA Q + dQ B - pdV 



(34) 



where U = J nu dV is the internal energy of the system, 
|42| , and as we can see, its change is due to a combination 
of heat Q added to the system and work done on the 
system. Furthermore, we have that 



dt 
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<j> n( V In n)' ■ n dS = <J> nv p ■ n dS, 



is the corresponding quantum heat flux due to the quan- 
tum nature of the KG equation. The second term on the 
right hand side of equation (IMl) would make the crucial 
difference between a classical and a quantum first law of 
thermodynamics. 

Analogously, for the magnetic heat we have 



dt 



' l<hl = [v-J B dV = jv (nj B ) dV 

= 7 J n (~W + (v ' V)A + (A ' v),y ) ' vdv 



(35) 



where the vector potential A fulfills the Maxwell equa- 
tions, in terms of the fluxes (ITB)) it reads 



piMu 



-3 



(36) 



where as usual F^ v — d p A v — d v A^. In terms of the 
vector and the electric potential, the Maxwell equations 
are given by 



UA = j 

2en 

Dip = -j -771 



(37a) 
(37b) 



where we have used the Lorentz gauge. Observe that 
the fluxes contain the information of the velocity of the 
fluid and of the electromagnetic term as well. This point 
will be important for the superconductivity. 



THE CONDENSATION TEMPERATURE 

First, we analyze the condensation temperature in the 
ideal case. In this situation the associated single-particle 
dispersion relation is given by 



E 2 = (pc) 2 + {mc 2 f 



(38) 



In order to compare our case with well-know critical con- 
densation temperatures, we start analyzing the ultra- 
relativistic and the non relativistic cases. For these we 
can express the dispersion relation (|3"8")l as E ~ p s [io| . 
where s — 1 corresponds to the ultra-relativistic system 
and s — 2 stands for the non-relativistic case (these top- 
ics have been extensively studied, see for example [34J439| 
and references therein). Taken into account the number 
of antiparticles, in the ultra-relativistic limit the conden- 
sation temperature is given by 



3h 3 c 



where 71 is the charge density 71 — 



1/2 



N—N 



(39) 



, iV is the num- 
ber of particles, N corresponds to the number of anti- 
particles and V is the volume of the system. The quantity 
N — N can be written as 



N - N 



^2[n p (fi,T) -n p (/i,T)] 



(40) 



where 71 obeys the Bose-Einstein distribution function 

1 



oP{E-v) _ 1 ' 



(41) 



Here fj, is the chemical potential and j3 — l/fc^T. Simi- 
lary for the anti-bosons n we have 



1 



(42) 



The chemical potential is bounded as \fi\ < rac 2 and at 
the condensation temperature we have that = mc 2 , 
which corresponds to the minimum of the associated en- 
ergy. For ultra-relativistic systems without antiparti- 
cles [4lj ]. the corresponding condensation temperature is 
given by 



knT r 



( 



ir 2 N \ 1/3 
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Finally, in the non-relativistic case 4(| 41 



\VC(3/2)J 



2irh 2 



m \V((3/2) 



2/3 



(43) 



(44) 



where £(x) denotes the Riemann zeta function. Expres- 
sions (f39|). (|43| and ((44)) corresponds to 3-dimensional 
systems, but these expressions can be generalized to dif- 
ferent dimensions. 
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On the other hand, the analysis of ideal and weakly 
interacting non-relativistic Bose-Einstein condensates 
with a finite number of particles, trapped in different 
potentials (see (il. [20I- 33. 4jJ-42| and references therein) 
shows that the main properties associated with the con- 
densate, in particular the condensation temperature, 
strongly depends on the characteristics of the trapping 
potential, the number of spatial dimensions and the as- 
sociated single-particle energy spectrum. 

In what follows we obtain the condensation temper- 
ature associated with our system within the semiclassi- 
cal approximation 0, [2(J 42 1. Inserting plane waves in 
the KG equation Q , neglecting the term proportional to 
T 4 in (||)]), assuming that the temperature is sufficiently 
small and considering the low velocities limit, allows us 
to obtain the single-particle dispersion relation between 
energy and momentum [l9j 



2m 

mc 2 (f> 



A 



2mc 2 
e 2 A 



A 



4mc 2 



2mc 2 



(k B T) 2 (45) 



$(f,t) 



We remain that in this work we interpret k 2 

as the spatial density n(r, t) of the cloud, being k the 
scale of the system. Notice that if we set A = 8irh 2 c 2 aK 2 
and <f> = ar 2 in (l45j) . where a = l/2(a>o/c) 2 and uj 
is a frequency, we obtain the semiclassical energy spec- 
trum in the Hartree-Fock approximation for a bosonic 
gas trapped in an isotropic harmonic oscillator [a 0, 42 1 , 



but with two extra terms due to the contributions of the 
thermal bath and to the electromagnetic field. 

In this work we consider only two cases, the first one 
is for A 2 proportional to a constant. As we will see, in 
this case the associated constant can be absorbed by the 
chemical potential u. For the second case we consider a 
dependence of the form A ~ r 2 . 

Assuming static thermal equilibrium n(r, t) ps n(r) 
[H, thus 



= K- 2 



n(r). 



(46) 



The spatial density within the semiclassical approxima- 
tion reads H [H 



n(r) = 



1 



d 3 pn(r,p), 



(47) 



(2tt/i) 3 

where n(f,p) is the Bose Einstein distribution function 
given by [8j, |20| 

1 



n(r,p) 



(48) 



The number of particles in the 3-dimensional space obeys 
the normalization condition [8l, L2Cj| 



N = I d 3 rn(r). 



(49) 



Integrating (14"T|) over the momentum space allows to ob- 
tain the spatial density associated to this system 



n(r) 



where Z — exp[/3(/i- 



mk B T 
2nh 2 



3/2 



53/2 (Z) 



(50) 



f4n(r1 - ^P^- ~mc 2 (f>~ f^)l. 

2mc z V I Amc z " 2mc 2 /J 

The function q v (z) is the so-called Bose-Einstein func- 
tion defined by [401 ] 



1 



r(iz) 



l dx 



1 



(51) 



being T(v) the Gamma function. In order to calculate the 
condensation temperature let us suppose that our gas is 
trapped in a harmonic oscillator type-potential <fi ~ r 2 . 
Clearly, this can be extended to a more general poten- 
tials. Expanding (j5"0")l at first order in the coupling con- 
stant A, using the properties of the Bose-Einstein func- 
tions [40(, allows us to express the spatial density as fol- 
lows 



n(r) fa n (r) ~ \g 3 / 2 {z{r)) 



2mc 2 n B T 



5l/2 (z(r))(52) 



i k b T gi/ 2 {z{r)) 
4mc 2 5 3/2 (z (r)) 



where 



no(r) = A 3 g 3 /2(z(r)), 



(53) 



is the density for the case A = 0, A = (2Trh 2 /thkbT) 1 ^ 2 is 
the thermal de Broglie wavelength, and z{r) = exp(/3(/i— 
amc 2 r 2 — e 2 A 2 /2mc 2 )). In the case A 2 ~ r 2 , with the 
help of the normalization condition (1491) we obtain 

3/2 



N 



) (k B T) 3 g 3 (e^) 
G 3/2 (e^) 



\2iih?) 

\ K - 2 m 2 (k B T) 7 / 2 



(54) 



167T 3 / 2 C 2 ftW/ 2 



A 



/ TO !/3 \ 3/2 



4:C 2 \2Qh 2 ) 



where 



G 3 /2(e Pfl ) ^2 jl/2 7 -3/2(j 1 j)3/2 ' 

i,i=l v ■' ' 



(55) 



and fl — amc 2 + const e 2 / 2mc 2 . From expression (|54l) . 
we easily obtain the case A 2 = const, 



N 



1 



2ac 2 ft 2 
Ak~ 2 



3/2 



(— 

V27T 1 / 



1/6 



(56) 



2C 5 V27T 1 /2^2 a l/2 
A / 1 \3/2 

Amc 5 \2ah 2 



{k B T) 7 ' 2 G V2 {e^-^) 



(k B T) 4 g 2 (e^-% 
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where 7 is defined as 7 = const e 2 /2mc 2 . We notice im- 
mediately from the expressions (|54l) and ((55)) that if the 
function A 2 is position dependent, the correction in the 
number of particles can be associated to an effective ex- 
ternal potential. On the other hand, when the function 
A 2 is constant, the correction can be associated to an 
effective chemical potential. If we set A = in expres- 
sions ([54")) and (1551) . we recover the expressions for the 
number of particles in the non-interacting case. In the 
thermodynamic limit, in the non-interacting case A = 0, 
at the condensation temperature the value of the chemi- 
cal potential is a = 8] . If we further assume that above 
the condensation temperature the number of particles in 
the ground state is negligible, this allows us to obtain an 
expression for the condensation temperature To. For the 
case A 2 ~ r 2 , we obtain 



N \i/3/2fi/i 2 \V2 



vC(3)y v m 

Analogously, for the case A 2 = 




(57) 



(58) 



Otherwise, using the properties of the Bose-Einstein 



-T/k B T 



when 



functions, particularly for (73 (e~ 
— 7/fcsTo — > [40(, where T is the condensation 
temperature associated to the case A 2 = const ^ 0, we 
obtain the shift in the condensation temperature respect 
to ([55)) caused by 7 7^ 0, in function of the number of 
particles, 



To -T Q 
To 



7 



C(2) 



3C(3) 2 / 3 W8o3 



(59) 



Clearly, the shift f|59[) vanishes when the number of par- 
ticles N —> 00 and tends to the value (|58p . In order to 
obtain the leading correction in the shift for the conden- 
sation temperature caused by the coupling constant A 
and the thermal bath in our system, let us expand the 
expressions ([54)) and ([56)) at hrst order in T = To, /i = 0, 
A = 0, and 7 = 0. Additionally, at the condensation 

temperature, the chemical potential in the semiclassical 

2 

?n(r = 0), as it is 



approximation is given by \x c 
suggested from expression ([5"0")) . thus 



2mc 



Ak- 2 to 1 /2( Kb t c ) 3 / 2 C(3/2) 



A 3/2 



2(2tt) 3 / 2 c 2 ?i 3 
V2nn~ 2 {K B T c y 
(2irc 2 h 2 ) 3 / 2 



(60) 



Expression ([SUf basically corresponds to the definition of 
the chemical potential at the condensation temperature 
in the usual case @, [2(| , except for the extra term contri- 
bution due to the thermal bath, and comes from the fact 
that g 3/2 (e~ s ) ~ C(3/2) - |r(-l/2)|<5 1 / 2 , when S 



[40l | . Using these facts, we finally obtain the shift in the 
condensation temperature caused by A and the thermal 
bath, in function of the number of particles in the case 
A 2 ~ r 2 



where 



_ 1 / C(3/2)C(2) , 
Ul ~3C(3) I, 2(2tt) 3 / 2 ~ G3 / 2 ^)' 



e 2 



1 (2A)V 2 C(2)7r 
3C(3) V4mc 2 (2tt) 3 / 2 k 2 /i 3 c 3 



and 5 = (2Qh 2 /m) 1/4 with T defined in l[57). A similar 
analysis, leads us to the shift in the condensation temper- 
ature caused by the coupling constant and the thermal 
bath associated with the case A 2 = const 



T c -Tq 
To 



at; 

To 



= -A 
+ 7 



,1/2 



K 2 h 3 C 2 

cm 



3C(3) 2 / 3 ficV8a 3 



6iEW 1/6 + \0 2 Z 2 N 1/3 (62) 
TV" 1 / 3 , 



with 5 = (2ac 2 fi 2 ) 1/4 and T defined in flSJJ). In the 
case A = we recover from ([^2^) the shift (|59l) . as ex- 
pected. From (I5T|) and (|6"2"|) we observe that the conden- 
sation temperature T c is corrected with respect to the 
usual case T as a consequence of the thermal bath and 
the contribution of the field A 2 , included in the semi- 
classical energy spectrum (|4"5j). Additionally, we notice 
that the critical temperature associated to the symme- 
try breaking (|10[) becomes very large when A — > and 
the condensation temperatures (|6"Tj) and (jUJ) tend to the 
non-interacting values (|ST)) and (|5"9l respectively. Set- 
ting a = l/2(wo/c) 2 and A = 8nh 2 c 2 K 2 a into expressions 
(|61[) and (|62[) we recover the condensation temperature 
for a bosonic gas trapped in an isotropic harmonic os- 
cillator potential, but corrected by the contributions of 
the thermal bath and the external field A 2 . For the sake 
of simplicity, let us analize the correction over the usual 
result caused by the thermal bath and the external field 
A 2 = const in (|6"2"|) . For instance, in the case of 87 Rb, 
with a - 10~ 9 to, N - 10 6 , and u ~ 10Hz, we ob- 
tain from the second right hand term in f[H2~[) a correction 
up to 7.9 x 10~ 78 k 2 + 7.5 x 10~ 38 k, and for the third 
right hand term, which is independient of the scale k, 
up to const x 10 2 . In other words, the scale k must be 
very large (up to 10 38 ) and the external field A 2 must 
be very weak (smaller than, or of the order of 10~ 12 ), at 
least near to the center of the system, in order to obtain 
relevant corrections over the usual result under typical 
conditions. For these values the symmetry breaking tem- 
perature is ~ lO 6 -^. Finally, with the experimental data 
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given above, we obtain for the first right hand side term 
of expression (|^|) the usual shift ~ 10 -2 , as expected 
[20l j . The temperatures T^ B and T c (or more specifi- 
cally, AT c /Tq) are related through the coupling constant 
A, this fact could be used as a criterium to compare both 
temperatures and in principle, to infer bounds related to 
the scale k. 



THE PHASE TRANSITION 

From hereafter we study the transition between the 
$ = state to the minimum R min = k B /2y/(Tj? B ) 2 - T 2 
with T < T? B . 

During the time when T >> T^ B there are not scalar 
particles in the symmetry broken state. Below the critical 
temperature T < T^ B , close to the local minimum the 
density oscillates around the value n = k^K 2 ((T^ B ) 2 — 
T 2 )/4. We study the case when the function S in (fT4|) 
has the simple expression, S = sot, with sq << mc/H in 
the non-relativistic limit. This implies that the velocity 
v = as well. If there does not exist an external force in 
the system, then = 0. In this case the viscosity (dissi- 
pative term) of the BEC might in fact contain the whole 
information of the phase transition. From equation (|23[) 
we observe that the viscosity Vcr contains only a therm 
with the anisotropics of the temperature. That means 
that when the temperature of the system isotropies the 
fluid becomes a superfluid. Furthermore, from the ex- 



pression for flux (|16al) we observe that the vectorial flux 
contains only a term with the vector potential A. Thus, 
the flux expression (|16a|) becomes the London equation, 

2ne 2 . 



indicating that the system becomes superconductor. The 
Maxwell equation (|37a[) becomes 



UA = 



2ne 2 



which is the Proca equation, indicating that the photon 
acquires a mass 2ne 2 /n 2 c. Obviously, playing with the 
conditions of the system we can find situations with su- 
perconductivity or superfluidity in different situations. 

Finally, to illustrate the previous exposition, we give 
the following example. Suppose that in the system there 
are only condensed and excited particles of the same 
specie. We interpret n = K 2 fc|((T ( f s ) 2 - T 2 )/4 as the 
particles density when the symmetry J/(l) has been bro- 
ken. At T = we expect that all particles have passed 
from the symmetry state into the broken symmetry one, 
such that the total particle density in the system is 
tit = K 2 fc| 3 (T ( f B ) 2 /4. Thus, in any time the number 
of particles Nq in symmetry broken state is given by 

2 

/ / 1 

1 - 



N Q =N 



( t sb) 



(63) 



where n/nr = Nq/N. Note that in this case the ex- 
ponent 2 in the critical temperature appears naturally. 
As usual this expression shows the dependence of the 
condensate fraction No/N as a smooth function of tem- 
perature from T ~ T C SB down into T = 0. In this case, 
the finite temperature terms are obtained from the one 
loop corrections of the SF density and are in complete 
agreement with the standard theory, [H, [H, [3^ ■ 

Observe that the total particle density jit = 
K 2 k%(T^ B ) 2 /A of the scalar particles reaching the sym- 
metry broken state at T = can only be determined 
experimentally and fix the value of the scale k. So in 
principle we are able to mimic the result that in the pres- 
ence of interactions we have Nq < N. The main idea we 
want to point out here is that these phenomena might be 
equivalent for a BEC on earth as well as for cosmological 
scales, so this model might be tested in the laboratory. 
If confirmed, the phase transition of a BEC can be ex- 
plained using quantum field theory in a straightforward 
way. 



CONCLUSIONS 

In this work we studied the phase transition of a boson 
gas with zero spin represented by the KG equation with 
a Lagrangian containing a U(l) symmetry, with mass 
m and self-interaction A, given by a mexican hat SF 
potential, immersed in a thermal bath at temperature 
T, close to the critical temperature of symmetry break- 
ing, up to one loop in perturbations theory. We rewrite 
the KG equation and interpret it as a generalized Gross- 
Pitaevskii one at finite temperature. We show that the 
transition from the phase with the U(l) symmetry to the 
phase with this symmetry broken can be interpreted as a 
phase transition from the gas state to the condensation 
state of the Bose gas. We obtained the condensation 
temperature as well. By rewriting the generalized Gross- 
Pitaevskii equation in terms of hydrodynamic quantities 
we were able to derive the thermodynamic of the phase 
transition, with viscosity and dissipative terms and find 
that the first law of the thermodynamics contains a new 
term that is a direct consequence of the quantum char- 
acter of KG equation. The main result of the present 
work is that the phase transition do not take place at the 
same temperature and conditions of the condensation. 
We saw that, for example, Bose-Einstein condensation 
take place in a 87 Rb crystal in normal density conditions, 
but that the phase transition can take place only in very 
hight density conditions, for example, at densities like in 
a neutron star. On the other hand, other materials have 
phase transitions and/or condensations in different con- 
ditions of density and temperature, depending on their 
mass, self-interaction parameter A and the value of the 
scale k. We gave different examples. 

It remains to see whether this generalization of the 
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Gross-Pitaevskii equation can describe the transition of 
a Bose gas into a BEC state in the laboratory. In other 
words, we propose that the superfluid and/or supercon- 
ductor like-behavior in a BEC can be measured experi- 
mentally in a laboratory, in order to compare the results 
given here with realistic systems. 

This work was partially supported by CONACyT 
Mexico under grants 166212, 132400 and 10101/131/07 
C-234/07 of the Institute Avanzado de Cosmologia (I AC) 
collaboration (http: / /www.iac.edu.mx/) . 



[1] 

[2] 

[3] 

[4] 
[5] 
[6] 

[8] 



[9] 
[10] 

[11] 



[12] 
[13] 
[14] 

[15] 



Electronic address: tmatos@fis.cinvestav.mx 
Electronic address: ecastellanos@fis.cinvestav.mx 
M. R. Anderson et al, Science 269, 198-201, (1995). 
Ph. W. Courteille, V.S. Bagnato, V.I. Yukalov, I. Laser 
Phys., 11, 659 (2001). 

M. B. Pinto, R. O. Ramos and J. E. Parreira, Phys. Rev. 
D 71 123519, (2005). 

L. Dolan and R. Jackiw, Phys. Rev. D 9, 3320 (1974). 
S. Weinberg, Phys. Rev. Lett. 9 3357 (1974). 
T. Matos and A. Suarez, Europhys. Lett. 96 (2011) 56005 
[arXiv:1110.3114 [gr-qc]]. 

T. Matos and A. Suarez, arXiv:1103.5731 [gr-qc]. 

C. J. Pethick and H. Smith, Bose-Einstein Condensation 

in Diluted Gases, Cambrige University Press, Cambrige 

(2006). 

E. Castellanos and T. Matos, arXiv: 1202.3793. 
Dany Page, James M. Lattimer, Madappa Prakash An- 
drew W. Steiner ApJ 155, 623, (2004) 
Juan Magana, Tonatiuh Matos, Victor Robles, Abril 
Suarez. A brief Review of the Scalar Field Dark 
Matter model. J.Phys.Conf.Ser. 378 (2012) 012012, 
arXiv:1201.6107 

O. Erken, P. Sikivie, H. Tam, and Q. Yang. Phys. Rev. 
Lett. 108, 061304 (2012). 

Juan Barranco Argelia Bernal. Phys. Rev. D 83, 043525 
(2011). 

Victor H. Robles, T. Matos. Exact Solution to Finite 
Temperature SFDM: Natural Cores without Feedback. 
arXiv: 1207.5858 

T. Chiueh, Phys. Rev. E 57, 4150 (1998). 



[16 
[17 

[is: 

[19 
[20 

[21 
[22 
[23 
[24 
[25 

[26 
[27 
[28 
[29 
[30 
[31 
[32 

[33 
[34 
[35 

[36 

[37 

[38 

[39 
[40 

[41 

[42 



D. Bohm, Phys. Rev., 85 (1952) 180. 

G. Grossing, Phys. Lett. A, 388 (2009) 811. 

Oliver, X. and de Saracibar, C. A., Mecanica de Medios 
Continuos, edited by Vieira, E. and Car E. (UPC) 2000. 
B. Zwiebach, Cambrige University Press, Cambrige 
(2004). 

F. Dalfovo, S. Giordini, L. Pitaevskii, S. Strangari, Re- 
views of Modern Physics, Vol. 71, No. 3, April (1999) pp. 
463-512. 

Ketterle, W., and N. J. van Druten, Phys. Rev. A 54 
(1996b). 

Bagnato, V., D. E. Pritchard, and D. Kleppner, Phys. 
Rev. A 35 (1987). 

S. Grossmann and M. Holthaus, Phys. Lett. A 208 
(1995). 

Giorgini, S., L. Pitaevskii, and S. Stringari, Phys. Rev. 
A 54 (1996) 

H. Haugerud, T. Haugset, F. Ravnal, Phys. Lett. A 225 
(1997). 

H. Shi and W. M. Zheng, Phys. Rev. A 56 1046, (1997) 

Z. Yan, Phys. A 298, 455 (2001) 

L. Salasnich, Int. J. Mod. Phys. B 16, 2185 (2002) 

O. Zobay, J. Phys. B 37, 2593 (2004). 

V. I. Yukalov, Laser Phys. Lett. 1, 435-461 (2004) 

V. I. Yukalov, Phys. Rev. A 72, 033608 (2005). 

A. Jaouadi, M. Telmini, and E. Charron, 

arXiv:1011.6477vl [cond-mat. quant-gas] (2010). 

V.I. Yukalov, Phys. Part. Nucl. 42 (2011) 460-513. 

F.Jutter, Z. Phys. 47, 542 (1928). 

H. E. Haber and H. A. Weldon, Phys. Rev. Lett. 46, 
(1981) 

H. E. Haber and H. A. Weldon, Phys. Rev. D 25, 502 
(1982). 

M. Grether, M. de Llano, and G. A. Baker, Phys. Rev. 
Lett. 99, 200406 (2007). 

S. Singh and P. N. Pandita, Phys. Rev. A 28, 1752 (1983) 

S. Singh and R. K. Pathria, Phys. Rev. A 30, 442 (1984). 
R. K. Phatria, Statistical Mechanics, Butterworth Heine- 
man, Oxford (1996). 

W. Greiner and Neise Stocker, Thermodynamics and Sta- 
tistical Mechanics, Springer- Verlag, New York (1995). 
Pitaevskii, L. Stringari, S. Bose-Einstein Condensation. 
Ed. Birman, J. et. al Clarendon Press, Oxford, (2003). 



